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Abstract 

We consider extensions of the standard model inspired by intersecting D-brane constructions, 
in order to address flavor mass textures. We include additional anomalous gauge symmetries, and 
scalar fields to break them and to generate Froggatt-Nielsen mass terms. Green-Schwarz axions 
are included to cancel mixed anomalies rendering the models consistent. At low energies, a resid- 
ual anomalous global symmetry remains, and its associated pseudo-Goldstone mode becomes the 
physical axion, which can be interpreted as an axion arising from open string modes. General 



OA 



considerations show that such axions are very common in D-brane models and can be completely 
incompatible with current bounds. Astrophysical constraints are placed on the axion both by 
including neutrino masses in the Froggatt-Nielsen scheme and considering QCD instanton contri- 
butions to the axion mass. We find simple models where the axion decay constant is in the allowed 
range, but only one such minimal model with this property is free from excessive fine tunings 
elsewhere. We also note that generically addressing flavor textures for the CKM matrix leads to 
deconstructed extra dimensions. 

PACS numbers: 11.25.Uv,14.80.Va,11.30.Hv 



I. INTRODUCTION 



The standard model includes no fundamental explanation for the hierarchy of elementary 
particle masses. The masses, and the Yukawa couplings which determine them, are spread 
across five orders of magnitude. This is the problem of flavor in the standard model of 
particle physics. Here we also include the fact that the CKM matrix is an input to the 
standard model and is not the result of a calculation. This problem of flavor begs for 
an explanation from a more fundamental theory of physics. A natural candidate for such 
a theory is string theory, which in principle would also unify the gauge interactions with 
gravity. As the standard model is weakly coupled at high energies in current accelerator 
data, it makes sense to take a weakly coupled string theory model to describe the standard 
model physics plus its extension to flavor physics. For the purposes of this paper we will 
work in models where the standard model is realized on stacks of D-branes and instead of 
a top-down approach, we will assume that everything can be described in terms of effective 
field theory of open string states of spins zero, one half and one, plus massless closed string 
states as required for anomaly cancelation. All these states are effectively massless at the 
string scale: we expect all masses to be generated at low energies by field theory mechanisms. 

A natural way to explain a pattern, such as the flavor problem, is in terms of symmetry 
and a symmetry breaking order parameter. This would explain the hierarchies by having 
different operators in the standard model transform under different representations of the 
symmetry group. If the symmetry is continuous, such a symmetry breaking would lead 
to a massless Goldstone boson as a low energy remnant. Such a possible remnant can be 
an axion and it is interesting to find models where the axion is related to this symmetry 
breaking of a flavor symmetry. It has been argued that within string theory and perhaps 
in any theory of gravity, there can be no global symmetries ( see jl| for a modern analysis 
of this argument), so any symmetry solution that is compatible with string theory should 
also come with associated gauge particles. In principle the symmetry breaking process 
would render them all massive and invisible at low energies. However, we will see that 
if the symmetry we gauge has mixed anomalies, there is a nice way out of this problem 
that can lead to a viable axion: the idea is that the symmetry can become an approximate 
global symmetry of the low energy theory that is only broken by non-perturbative effects. 
This was first found by Ibanez and Quevedo |J] in D-brane models and they argued this 



way that the proton would be long-hved. This survival of the gauge symmetry as a low 
energy global symmetry remnant in perturbative effective field theory makes it possible to 
generate approximate symmetries like those that determine the axion dynamics and hence 
has observable consequences. This is what we will explore in this paper, in the context of 
D-brane inspired models. It is important to notice that even though we will be working with 
string inspired models, the arguments we give can be applied in more general setups and can 
be argued completely within effective field theory. The string putative embedding restricts 
the choices of examples we work with giving a simplified set of rules to build models, as 
well as specific predictions within the models that can be compared to experimental and 
observational bounds. Many of the models we explore will be ruled out this way, but not 
all of them. More general effective field theory models built with similar rules will have 
many more free parameters related to mixings that are forbidden in brane models, as well 
as charges not being quantized by simple rules and are much less predictive. 

Let us begin our discussion with having a flavor symmetry that is gauged to address the 
problem of flavor. This is not a new idea. The Froggatt-Nielsen mechanism provides 
a way to generate the hierarchy of masses while maintaining 0{1) Yukawa couplings, by 
adding a gauged flavor symmetry to forbid some tree-level fermion masses. This symmetry 
is spontaneously broken, and effective nonrenormalizable operators of various dimensions 
provide fermion mass terms. These operators contain different powers of the order parameter 
for symmetry breaking, which in turns generates hierarchies between fermion masses. 

In this paper we make use of this mechanism in the context of simple string-inspired 
models. This paper is a natural extension of the models introduced in j4|. The abelian 
Froggatt-Nielsen gauge symmetries in these models can have mixed anomalies, which are 
then cancelled by a generalized Green-Schwarz mechanism jsl. This requires the introduction 
of a massless pseudoscalar boson that couples to the anomaly and cancels it. Thus, these 
particles couple to the Standard model like axions and we will call them axions. When these 
axions cure the anomalies they end up providing Stiickelberg mass terms for the anomalous 
gauge bosons. Thus they are eaten up and disappear from the low energy effective theory. 
However, it turns out that in D-brane models the global parts of these gauge symmetries 
survive in the low-energy theory ; we will refer to this phenomenon as the Ibanez-Quevedo 
mechanism. It should be noted that typically, all symmetries in string models should be 
gauged, and finding systems with approximate global symmetries usually requires unnatural 



fine tuning; the Ibaiiez-Quevedo mechanism provides a neat way to avoid this. In order to 
generate fermion mass terms in our models, these residual approximate global symmetries 
from the Ibafiez-Quevedo mechanism must be broken. We identify these residual global 
symmetries as Peccei-Quinn symmetries. 

The Peccei-Quinn symmetry is a classic mechanism which addresses the smallness of the 
QCD vacuum aug^ fl. Oue of its low energy cousequeuees is an axion whose properties 
can be predicted [7|, |8[. The simplest types of axion have been ruled out, meaning that 
the Peccei-Quinn symmetry has to be broken at scales much higher than the electroweak 
scale, producing an invisible axion. Experimentally, there is a narrow 'axion window' which 
constrains the decay constants of these particles to lie between 10^ and 10^^ GeV. These 
lower and upper bounds on the scale of Peccei-Quinn symmetry breaking arise respectively 
from astrophysical and cosmo logical considerations; this is reviewed in 9|. 

Axions are a generic feature of stringy models, usually arising from closed string modes 
in the context of the Green-Schwarz mechanism. Typically, their decay constants are close 
to the string scale. Finding realistic axions in stringy models is not always easy, but one 
can not rule them out either flo|. Our models include axions of this type, but we include 
scalars arising from open string modes to break the global Peccei-Quinn symmetry, which 
induces additional axions. It is these 'open string axions' which survive as the physical 
axion at the lowest energies in our models. This feature has been found in some top down 
compactifications leading to supersymmetric models previously ll| (these models have been 



analyzed further in 



12|). 



The precise string inspiration for this model comes from studying intersecting D-brane 



setups (see 



13| for a review). The matter content of these models is usually presented as a 



quiver (or moose) diagram. We will follow these conventions. To simplify the analysis we will 
insist on requiring only the minimal particle content that is able to accommodate the flavor 
hierarchy generation mechanism at the level of effective field theory. The model we write is 
closely related to the 'Minimal Quiver Standard Model' introduced in pL4i]. The effective field 
theory will have high dimension operators suppressed by some high scale M. For simplicity 
we will call such scale the string scale, although the model does not need to have a precise 
string theory origin. We find that if the neutrino masses are included in the analysis, we are 
able to determine both M and the Froggatt-Nielsen (Peccei-Quinn) symmetry breaking scale, 
and this gives a prediction of the energy scale of the axion decay constant. Furthermore, we 



find that in one of the models the axion decay constant falls exactly in the allowed window, 
while other models can be ruled out because the corresponding axions are not allowed or 
there are some other basic problems with their phenomenology. 

Related work along these lines has been done. In [isl] a similar scheme of incorporating 
axions into a Froggatt-Nielsen mechanism is discussed. The authors use a rather different 
device to break the PQ symmetry in a supersymmetric setting, resulting in Higgs-axion 
mixing and axion masses and couplings which are not strongly determined by the axion 
decay constant. A similar result occurs in the models studied in Q. These features 
are not exhibited by the physical axion appearing in our models. Also, our gauged flavor 
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lich gauge nonabelian flavor symmetries can also yield 



symmetries are abelian, but mode 
axionic degrees of freedom, as in 

II. MINIMAL MODELS 

As we discussed in the introduction, we want a D-brane inspired model that contains 
the standard model and some gauged flavor symmetries. Considerations of orientifolded 
D-brane embeddings of the theory require that all SU (N) gauge symmetries be enhanced to 
U (N) symmetries at the string scale, but also make the inclusion of Sp{N) gauge symmetries 
natural. Thus, we require a U{3) color symmetry, but are also allowed to use Sp{l) ~ SU{2) 
rather than enhancing isospin symmetry to U{2), avoiding a possible extra U{1) gauge 
symmetry. Having a realistic model with minimal matter content requires an extra brane 
for a gauged U{1) symmetry, so that the theory has at least a f/(l) x f/(l) gauge symmetry. 
This is true regardless of whether we have SU (2) or U (2) associated to the weak interactions. 
Although this could in principle be avoided by using models that arise from a grand unifled 
theory, or that contain quarks in a two index antisymmetric representation of f/(3), they 
are extremely problematic as many important couplings of the standard model would need 
to be generated non-perturbatively ( see [3] for a discussion about this issue). 

The minimal matter content would have us select an SU{2) symmetry. As is typical 
of such setups, only hypercharge is a non-anomalous gauged symmetry. The other U{1) 
symmetry has mixed anomalies that are cancelled by the four dimensional version of the 
Green-Schwarz mechanism [Sj]. This produces a heavy Z', and it is possible to generate all 
couplings in the standard model at tree level, so that the only new physics is the presence 



of a heavy Z' . These considerations produce the minimal quiver model 14 1. 

Now, we want to make the model richer by adding a gauged flavor symmetry that dis- 
tinguishes various quarks from each other. The simplest such extension requires only one 
more f/(l) gauge field. One way to include an extra gauge field is to add another brane to 
the minimal quiver model, giving a four stack model rather than a three stack model (here 
stacks also refer to nodes in the quiver diagram). The gauge anomaly cancellation condition 
tr [T"{T^, T'^Y\ — can be used to fit the right handed fermions of the standard model with 
two up-type quarks and one down-type quark at one node, and the opposite arrangement at 
the other node. The second possibility which we will study further below is to use a f/(2) 
gauge theory for the weak interactions, rather than SU{2). This second choice produces a 
three stack model instead. 

At this stage, we need to fit the Higgs doublet. There are two choices we can make: we 
can have a model with a Higgs field associated to one U{1) node but not the other; either 
choice forbids Yukawa couplings for three of the quarks, and possibly all of the leptons. The 
three quarks for which Yukawa couplings are allowed are chosen to be the heaviest quarks 
with the correct quantum numbers. This is, either t, c, b or t, b, s. We will call these two 
choices model A and model B. This is where we need to add more fields to generate masses 
for the other generations. To communicate the electroweak symmetry breaking to the rest 
of the fermions that connect to the node where the Higgs field is absent, we introduce a 
new complex scalar field between the two U{1) nodes that is neutral under the standard 
model. Another possibility would ask us to have a two Higgs field model instead, one 
for each U{1) stack: such models besides having more particles and not being minimal 
would also typically lead to unacceptably large flavor changing neutral currents. Under our 
assumptions of minimality, we choose to have only a single complex scalar field joining the 
two U{1) stacks. We will call the complex scalar the Froggatt-Nielsen (FN) scalar. 
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TABLE I: Charges under the Froggatt Nielsen symmetry in the four-stack models. For the Higgs, 
we have two choices: we call the two models A, B 



FIG. 1: Quiver with extra scalar (j). The two choices for the Higgs doublet are labeled by the 
dotted lines and the model they represent h^'^ (for the models A,B, respectively). We have added 
arrows to indicate the U{l)ab charges of (/>. 

In the full model a U{1)b x f/(l)a x U{\)h is gauged. Hypercharge f/(l)y is the only 
anomaly-free linear combination of these. Only the right handed quarks and leptons are 
charged under ?7(l)a,6, which is the set of gauge groups under which the scalar is charged. 
Charges under these U{l)s are depicted in table [B Figure [1] shows the matter content in 
quiver form, with both choices for the (single) Higgs doublet present. In figured] we have 
left the edges corresponding to right-handed quarks unlabeled since their interpretation as 
light or heavy depends on which Higgsing scheme is used; generically, up-type right-handed 
quarks are represented by edges, while down-type quarks are represented by edges 
in the figure. Gauge invariant Yukawa or FN couplings are constucted by tracing closed 
loops in the diagram. Some of these are high dimension operators. The neutrino masses are 
also generated in this way, and depending on the choices that are made, they give rise to 
dimension 5 or dimension 7 operators if they require the introduction of FN fields to render 
the couplings that generate masses gauge invariant. 

A mass hierarchy from a FN symmetry can also be introduced in a three-stack model, 
once again by a slight modification of the minimal quiver model. This is done by replacing 
Sp{l) by U{2)yj^ which introduces another gauged U{1) boson. Tadpole cancellation requires 
that the three left-handed quark doublets be split into two sets, with two generations in the 
fundamental and one in the antifundamental representation of f/(2)u,. This also ensures 
that the FN symmetry acts differently on different families as we want for our models. 
Incidentally, this prevents the U{1) C U{2) from giving contributions to the hypercharge 



Y. Fixing the hypercharge of the Higgs doublets requires fixing the arrows of the field 
at node U{l)b- The choice now becomes on how the Higgs transforms under U{2)y^: in 
the fundamental or the antifundamental representation of U{2)y^. As far as the SU(2)w is 
concerned, both of these choices are equivalent, but they differ in how the Higgs is charged 
under the U{l)w. This is pictured in figure O As in the previous models, there are two 
ways to charge the Higgs doublet, resulting in two models which we refer to as model C and 
model D. 

Once again, some Yukawa couplings are forbidden, and we include an additional scalar 
in order to generate FN terms and fully communicate electroweak symmetry breaking. This 
scalar is chosen to be in the antisymmetric representation of U (2) so that is a singlet under 
standard model isospin, but remains charged under the extra U{l)w C f/(2)^. Left-handed 
neutrino masses are generated in identical fashions to models A and B: by closing loops in 
the diagram and paying attention that the arrows at the nodes U{l)h and U (2)^ are oriented 
so that the gauge charges cancel. 
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FIG. 2: A three-stack model yielding FN terms. Both choices for Higgs charges, models C and D, 
are depicted in the diagram. 
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TABLE IL Charges under the Froggatt Nielsen symmetry in the three-stack models. For the Higgs, 
we list the two choices C, D 



Given the quantum numbers we have described, the minimal set of couphngs that generate 
all Yukawa couplings, including masses for left-handed neutrinos, is shown in tables IllltlVIl 
The high dimension operators (we will call these the Froggatt-Nielsen terms of the action) are 
suppressed by some high scale, M, which we take to be the string scale; this is a minimality 
assumption in order to be predictive. The suppression splits the fermions into light and 
heavy tiers. The heaviest field in each of the light and heavy tiers is taken to have an order 
one dimensionless coupling to the scalar fields; the ratio between tiers is set by the mass 
ratio of the heaviest quarks in each tier. This can be used to get an order of magnitude 
estimate on both (0) and the string scale M by examining mass ratios. We assume that the 
heaviest neutrino mass saturates the upper bound of rriy^ ~ 2 eV, and using the usual Higgs 
VEV {h) = 246 GeV we get the estimates given in the tables in section |Vl 

Notice that with this procedure we are able to predict all relevant scales for the model 
from known data. The reason for this is that we have used the neutrino masses to find 
the scale of new physics, since these are the parameters in the standard model that have 
been observed that are sensitive to the high scale physics (it arises from non-renormalizable 
couplings at the electroweak scale in the absence of right handed neutrinos, or when these 
have been integrated out). If in a given model the neutrino masses have a different origin 
than the string scale or the Froggatt-Nielsen scale, the above constraints change and tend to 
push the string scale higher. At this stage we must ask what is the physical implication of 
these deductions for low energy physics. In particular, we must determine how (0) is related 
to the decay constant of the axion appearing in the low energy theory, since the stringent 
bounds on this decay constant can be used to rule out models of axions. In order to proceed, 
we need to examine the effective field theory more carefully. We will do this in two steps. 
In the next section, section IIIIl we describe how to extract the axion in the general class of 
models we have constructed. In section IIVI we determine the axion model parameters that 
are needed to understand the axion couplings to matter. 

Notice that our models generate various textures for quark masses. These have been 
analyzed in a fairly general way in [20] for other models (see also 21[ for specific examples), 
but the general form of the textures is very similar. 



III. THE PHYSICAL AXION 



Let us begin our discussion with the U{1) gauge fields. We have a kinetic term of the 
form 

and there are no mixed terms between the kinetic terms in D-brane models at tree level. 
These can be generated by radiative corrections, but we will ignore this at this stage. 

The U{l)i charges Qi have mixed U{lY, U{l)Sp{lY , and U{1)SU{?>Y anomalies in each 
of the models. The full anomaly needs to be cancelled in order to have sensible physics. In 
each of the four models, we find that one linear combination Qy of these charges is anomaly- 
free, and two are anomalous. In all four models, Qy couples to the anomaly-free gauge field 
y^, which is some linear combination of the original three U{1) fields. We take the two 
remaining, anomalous charges Qi to be orthogonal to the anomaly-free charge, and to each 
other. These linear combinations can be determined by the vanishing of the three-way f/(l) 
anomaly. Knowing all three anomalous and anomaly-free linear combinations allows one to 
determine the mixing parameters of these gauge fields (analogous to the weak mixing angle 
sin^ 9w) in terms of the three original couplings gi. 

Many mixed anomalies remain uncancelled. We have taken care of the QyQiQ2 anomaly 
already, by the precise choice of Q1Q2', the rest of the triangle diagrams involving the 
gauge bosons associated with Qi are cancelled via the Green-Schwarz mechanism. The 
remaining mixed U(l) anomalies can be cancelled by including Chern-Simons terms in the 
effective action (see [22] for a discussion). When describing the Green-Schwarz mechanism, 
the mixed abelian-nonabelian anomalies are cancelled by the introduction of two additional 
pseudoscalar fields with axion-like couplings to the standard model gauge fields. These 
pseudoscalar fields transform non-linearly under the anomalous U{1) gauge transformations. 
Their gauge invariant kinetic term is of the form 

^K,,{d,s^ - Z;){d,s^ - Zf,) (2) 

and the scale of K is close to M^. 

Since the scalars transform non-linearly, we can choose the unitary gauge where the 
are constant. In this gauge the are clearly seen to acquire a Stiickelberg mass of order 
gM , where g is of order ga,gb,9c which are small at high energies (the factor of g is the 



usual open-string-closed string mixing), so they can be integrated out of the low energy 
effective theory. When integrating the Z' fields, this process produces only current- current 
interactions for the f/(l) gauge fields that have been integrated out, but these currents 
are invariant under the anomalous global f/(l)s. Starting from this simple action, the low 
energy effective theory seems to have a global symmetry in the lagrangian, but in this 

case this symmetry is anomalous. This is a particular realization of the Ibanez-Quevedo 
mechanism This mechanism is usually used to guarantee a long lifetime for the proton, 
but it serves generically to produce approximate global symmetries in the infrared effective 
theory. 

One of these symmetries is baryon number (this is used to prevent proton decay), while 
the other symmetry is the Froggatt-Nielsen symmetry of our lagrangian. The field is 
then a scalar field that breaks the Froggatt-Nielsen symmetry, and we get an additional 
Goldstone boson for this symmetry breaking. Notice that if the \J{V)p^ gauge symmetry 
were non anomalous, this Goldstone boson would have been eaten by the Froggatt-Nielsen 
gauge field and there would be no low energy remnant. This is because having a Stiickelberg 
mass for the gauge field would not be required by anomaly cancellation. Also, since the 
symmetry that is broken is anomalous in the low energy effective theory, the symmetry is 
not a true symmetry of the lagrangian anyhow and non-perturbative effects will lift the 
symmetry producing a potential for the would-be Goldstone boson. Such terms can be 
traced to the original lagrangian and they require corrections to the potential of the form 
y(0exp(zs), 0* exp(— is)), as described in [l5|. These terms are not allowed by perturbative 
computations in string theory so they can only be generated by D-brane instantons (see 16 1 



for a review of instanton computations). This means that the would be Goldstone boson 



is very light and gives rise to a light invisible axion [23|. If we take the parametrization 
= rexp(i^), then Q is the phase of that has been identified with an axion field at 
intermediate energy. 

A more careful analysis of the lagrangian shows that the field Q will mix with the 
fields and that the vev of will also contribute to the mass term of the Z fields. If the Z 
fields have a bare mass generated at the string scale, the new contribution would be a small 
correction at the Froggatt-Nielsen scale, so the correction from mixing of scalars is small 
and Q plays the role of the intermediate energy axion. This is what we expect in the limit 
of a high string scale - to a good approximation the physical axion is just the phase of 0, 



and the axion decay constant is simply /„ = (0). 

This can be seen exphcitly by considering the Stiickelberg terms for both the Green- 
Schwarz axions and for 9, which cause these three scalars to mix. Rewriting the Stiickelberg 

terms for all the scalars in terms of Y^-, we find that disappears, since is not 
charged under Qy in any of the models. Altogether we have 



We see that the Goldstone boson mixes with the Green-Schwarz axions to form new 
scalars 

^ fsiQzi i , {(t>)Pzi a _ s' + e0 



A*zi \/l + 



for mixing parameters and couphngs P^i and ji^i = V^(7si^^0^~+T(^)^^0^ ~ M-\/T+~?. In 
the approximations we take P^i ~ g^i ~ C'(l), /gi ~ M, and e = {(f))/M <^ 1. These axions 
ai are eaten by the heavy gauge bosons Z^. The heavy gauge bosons acquire the masses 
Mi ~ from eating these scalars. The combination oi s^,6 orthogonal to ai is the physical 
axion a, surviving in the low-energy theory: 



a — 



— e[s + s ) 



So indeed, up to corrections of order {(f)) /M, the physical axion is just the phase of (f). 
Since the Higgs is also charged under the Proggatt-Nielsen symmetry, when the electroweak 
theory is broken, there is in principle a further mixing between and the neutral would-be 
Goldstone boson from the Higgs doublet. However, at the level of the effective action at the 
weak scale no such terms seem to appear in the effective potential in our models. Any such 
mixing would be produced from mixing with the high scale Z' and is very suppressed, by 
v'^/M'^, so it can be safely ignored. 

The fields are as in the standard model, but Z differs somewhat. In models A and 
B, the situation is as follows; that for models C and D is analogous. Representing the Higgs 
doublet as h = {{h) + H + iz, w^), we find the Stiickelberg term 

(d^z + gj¥l{h) + 2g,{h){P,,.Z\ + P,,.Z\ + P^Y^)f 



upon writing in the YZ^''^ basis. We can determine the mixing parameter P^y by com- 
paring the abehan parts of the gauge covariant derivatives in the A"''^''^ and YZ^ bases, using 
the known relation between U{1) charges Qy = |(<5c — SQa — SQ;,), and the fact that the 
mixing matrix Pij is orthogonal. We find that Qb^li = ~^9yY^ + . . ., so that the relevant 
part of the covariant derivative is 



g^W^^T^ + g,AlQ, = {g^W^ - gyY,) + 2gb{Pb,iZl + A.^Zj) 



= MzZ^ + M,.Zl + M,2Zl 

The first term in parentheses in the center is the usual Z boson, and the other two terms 
form some combination of bosons. None of the Stiickelberg terms left after EWSB involve 
any mixing between a and any of the Goldstone modes from the Higgs (but there is mixing 
of these with the other 'axions' ai). Thus, this model does not provide an 'axi-Higgs,' as in 



e.g. 



15|- 



Since axions are so light, and can in principle be produced in profusion in low-energy 
nuclear processes, they must not couple strongly enough to nuclear matter to quench stellar 
fusion. This is the source of the lower constraint on the axion decay constant. This bound 
is mo del- dependent, so it is useful to check that the nominal axion window holds for the 
models under consideration here. In particular we wish to determine whether there are low- 
energy couplings which accidentally lead to fine-tuning and radically move the lower sill of 
the axion window in these models. 



IV. AXION COUPLINGS 

In order to determine the couplings of the axion to matter at low energies, we follow the 



procedure of [28[, which involves a series of matter field redefinitions. We perform two phase 
rotations of the matter fields, one to realize the PQ symmetry nonlinearly for physics below 
the PQ scale, and the other to eliminate axion-gluon couplings for physics below the QCD 
scale. These field redefinitions result in derivative couplings of the axion to currents, as is 
characteristic of Goldstone bosons. These phase rotations are also anomalous, resulting in 
axion-field strength couplings in the effective lagrangian. At the very lowest energies, we 
use these couplings to place bounds on the axion decay constant from nuclear processes. 



The first field redefinition is performed below the PQ scale, where in our models, the 
extra U (1) gauge symmetries are broken, and only the standard model gauge group remains. 
In this field redefinition, we rotate the phase of onto the fermion fields in the FN terms, 
making use of the unbroken PQ symmetry: 

(t)ilnph (e'"^^/-^0) (e'"^*/-^V) (e'^^^^-^V') h 

where -0, ip are generic left-handed Weyl fields; the bar indicates an antiparticle. We require 
the PQ charges Xi to sum to —1 in each FN term in order to cancel the phase of 0. For 
now, we leave these charges unspecified; we will discuss the choice of these charges below. 
This results in the coupling of the axion to the PQ current 



/a ' V 2 2 
The vector part of this coupling violates CP - we will discuss a resolution below. Since 

the PQ transformation is anomalous, it induces the axion-gauge interaction terms 



2 V 8^Va 

glNGG + glCwWW + glCyYv) 



167r2/a 

where the sum runs over the standard model gauge groups, and the trace runs over the 
fermion fields coupling to each. The extra overall factor of 1/2 arises because the field 
rotations are chiral. We have defined the PQ color anomaly = \ii[QpQ\. Below the 
weak scale the hypercharge and weak anomalies merge to give the electromagnetic anomaly 
coefficient E = Cw + Cy = ^"^^IQpqQeQe]-! so that only axion couplings to gluons and 
photons remain. 

The second field redefinition is performed below the weak scale, in order to eliminate the 
axion-gluon coupling in favor of couplings of axions to light quarks. In a sense this is the 
reverse of the first field redefinition, except that here we only rotate the fight quark fields 

The charge Qr — Qv + QaIs is anomalous, so its axial part modifies the effective axion- 
gauge interactions to 



QiV - 2N^tr[QA]^ aC^G' - {^E - 2iVtr [QaQeQe]) aFF 



with Qe the electric charge matrix of the hght quarks. We need trC^A = 1/2 in order 
to cancel away the axion-gluon coupling, but the choice of Qr is otherwise unrestricted. 
Convenient choices can be made to simplify various interactions. The light quark interactions 
are also modified to 

d,,a _ „ f /X^ - X;^ \ /X^ + 



The axial part of this current contributes to the axion-hadron couplings. The axion 



baryon couplings are determined via current algebra in 29| to be 



2(trQ^T") (Ftr57'^75[T^5] + Dtr57^75{T^5}) + -{tTQA)tTBr^,B 



fa { ' ^ , , 3 

where B is the baryon octet matrix, Q.A = ~ ^ ~ 0-44, D = 0.81, S 



0.13 are baryon decay parameters 30, l3l|]- This gives a derivative pseudovector coupling 
9aNN{dt,a/ fa)NYi-,N to nucleons. 

It remains to specify the charges Qpq,Qr for each model. One might hope that this 
could be done in such a way as to cancel off the vector part of the axion current, since it 
seems to be a source of additional CP violation in the strong sector, and complicates the 
calculation of bounds on the axion decay constant. This is not possible in general, because 
of the restrictions on the choice of the PQ charges. In these models the PQ symmetry is 
the global remnant of a gauge symmetry, so it must be some linear combination of the three 
U{1) generators in the quiver; that is, X^ = J2 O'iQiijP)- Given the f/(l) charges of the light 
quarks, and the fact that Qa must trace to 1/2, we find that it is not possible in any of the 
models to find a choice of Oj and Qv leading to (X^ — X:^)/2 + Qy = for each light quark. 

The vector part of the axion current can, however, be redefined away by yet another chiral 
rotation of the quark fields, done independently on each light quark. This has the effect of 
absorbing the offending CP violation into CKM matrix elements, so the seeming strong CP 
violation is really a weak CP violation. This might seem to violate electroweak symmetry, 
since up and down type quarks in the same isospin doublet are transformed differently, but 
this transformation is done in an effective theory at a scale where the electroweak symmetry 
is already broken. 



Given this, we choose Qpg and to simphfy calculations, and determine the couplings 
of the axion to matter in each model. Some of the Yukawa couplings must be somewhat 
finely tuned in order to reproduce the observed mass hierarchy. is chosen to simplify 
hadronic axion interactions as much as possible. In general, axions suffer mixing with tt^,"!]. 
Mass mixing can be eliminated by choosing 



2trM- 



2{1 + z + w) 



( 



\ 



\ 



w 



with Mq the light quark mass matrix and z = 171^/171^, w = rriu/ms. This choice of Qa can be 
seen as either eliminating axion-meson bilinears in an effective meson lagrangian including 
axions or directly via current algebra as eliminating off-diagonal terms in the Dashen 



formula for meson masses 



27 



331]. The resulting axion mass is 



^mif^z 



il + z)il + z + w)f^ 

If we also make the choice Qa = Qv then is purely right chiral, which simplifies the 
axion weak interactions. The Higgs is uncharged under the PQ symmetry, thus preventing 
axion-Higgs mixing. There is still kinetic axion-meson mixing, but these cross terms are 
suppressed by f^/fa- For each model, Qpq is chosen so that tr[(5pQ] = 1. In each case, 
the quiver structure allow us to choose Qpq to be a very simple linear combination of the 
anomalous U (1) charges; the specific choices for each model are described in the next section. 

We can now determine the axion-photon couplings. With the above choice of Qa, the 
axion-photon interaction term becomes 



E 

N 



24 + z + w 



aFF = gay aE ■ B 



167r2/a 31 + Z + W^ 

In all our models = 1/2, and E will be listed in the next section. 



V. MODEL FEATURES 



A. Four-Stack Models 

1. Model A 

Model A has the most natural set of FN couphngs for quarks; since the u, d, s quarks 
have suppressed masses in this model, the FN couplings must be down by only 10""^ at 
worst for the up and charm quarks. This does not work for the leptons, which naturally 
acquire masses near the top quark mass in this model. We consider the possibility of fine 
tuning at this stage (assuming that there is some other effect making these couplings small 
that has not been included in the model because of our minimality assumptions) and will 
check the axion phenomenology to see the features that show up. For this model, we choose 
QpQ — —Qb- 

The linear combination Qy — | {Qc — ^Qa — ^Qb) is anomaly-free, and corresponds to 
the electroweak hypercharge. This anomaly-free charge must be orthogonal to the two 
anomalous charges, which constrains them to have the form Qi = {'^i/^ — Pi)Qa+ PiQb+liQc- 
We wish to determine /3i,7i. First, we require the three-way U{1) anomaly coefficient to 
vanish; for model A, this is tr[(5y(5i(52] = y7i72 — 6^i^2- We also demand that the 
anomalous charges be orthogonal to each other. The overall normahzation of the anomalous 
charges is arbitrary, so finally we may choose — 1. This gives us two equations for the 
two unknowns /5j 

tr[gygiQ2] = y - = 

Qi-Q2^ {^-^^ Q-/32)+/3i/32 + l-0 
which gives us the anomalous charges 

Qi — Qa Qh + Q2 = ; Qa — Qb-\ — ^ ; ^-Qc 

6 6 5 5 

The following fermions receive their masses from dimension 4 operators: b, c, i. To 
accommodate the neutrino masses we find that we need dimension seven operators in the 



Dimension 4 couplings 


h'^quj + hqD + /i/e 


High dimension couplings 


(t>h^qU + ^'^hqdi + {(phHf 




10^ GeV 


Af ~ — ^ 


10^ GeV 


9app 


0.92 


Sarin 


-0.58 


E 


-13/12 


Sa-y 


-5 X 10"^ GeV~i 



TABLE III: Model A 



effective field theory. This is due to the combinations of the FN charge of the lepton doublets 
and the Higgs. The results are summarized in the table IIIII 

The most important thing to notice is that the since the neutrino masses come from 
dimension seven operators means that the neutrino masses require the high energy scale to 
be low From table IIIII we see that none of the axion couplings are unusually small, so the 
axion window is not shifted from its nominal range. The axion in model A is ruled out, 
since its decay constant falls many orders of magnitude below the axion window. 

2. Model B 

Model B fares much better than model A. It predicts light m, rf, c quarks and light leptons; 
this puts the smallest couplings in each tier down by 10~^. Only the Yukawa coupling of 
the strange quark is unnaturally small. Since this model differs only in its Higgs sector from 
model A, the anomaly-free hypercharge and anomalous U{1) charges are identical to those 
in model A. For this model, we choose Qpq = Qa- 

Again, none of the couplings to matter turn out to be unusually small. Model B produces 
a harmless axion, since its decay constant falls squarely in the center of the allowed window, 
as indicated in table IIVI Model B would have been preferred anyhow because in that model 
the leptons have masses comparable to light quarks without fine tuning. 

Notice that as seen in table IIVI the neutrino masses arise from dimension 5 operators. 
Thus they are only suppressed by the ratio {h)/M relative to the top quark, and this pushes 
the scale M higher. Also, all leptons have their masses from dimension five operators so 



Dimension 4 couplings 


h^qtj + hqdi 


High dimension couplings 


(p^h^qui + (phqD + cphle + (/it/)2 


(6) ~ ™= 


10^^ GeV 


M ~ ^ 


10^3 GeV 


9app 


-0.440 


9ann 


0.58 


E 


-1/3 




-3 X 10^1^ GeV-i 



TABLE IV: Model B 



they are naturally light. The string scale is compatible with expectations from D-brane 
intermediate scale models. At this stage we can declare success for this style of model 
building. 



B. Three-Stack Models 



1. Model C 



This model is problematic: the FN terms split the quark isospin doublets between the 
light and heavy tiers, with ^ controlling the difference between the mass tiers and fails to 
explain why some of the quarks are light by much larger fine tunings than in our previous 
examples. Just based on this one should strongly disfavor models like this one. Despite the 
schizophrenic pattern of Yukawa couplings, which is no better than those of the standard 
model, model C does however have the advantage of suppressing lepton masses and sharing 
important characteristics with model B. 

The anomaly-free hypercharge is given by Qy = ^{Qc — ^^Qb)- The orthogonal anomalous 
charges are 

g _j_ '^Y^ 1 8 T\X2 1 

Qi = — Q — Qw — '^^^^ '^c, Q2 = — Q — Qw + ^^^5^ 

We choose Qpq = \{Qw+Qb) in this model. Model D has identical anomalous, anomaly-free, 
and PQ charges. 

Once again, table |V] indicates that the axion has ordinary couplings to matter, so we 



Dimension 4 couplings 


hQu + h^qjd 


High dimension couplings 


(phqiu + (l)h^Qd + (phHe + {hlY 




10^2 GeV 




10^3 GeV 


9app 


-0.192 


9ann 


0.187 


E 


55/24 


9a'y 


3 X 10"^^ GeV-i 



TABLE V: Model C 



conclude that the axion decay constant in model C is close to the upper sill of the axion 
window, but not greatly outside it, rendering it mostly harmless. Indeed, the upper sill 
of the axion window depends on typicality arguments of the axion field expectation values 
in the early universe and should not be used to discard models; it does require additional 
mechanisms in the early universe that prevent the axion from acquiring random expectation 
values when the axion field decouples from the rest of matter and overdosing the universe. 

2. Model D 

This model has even more problematic mass hierarchies than model C. After all, the 
lepton masses are not suppressed and we are required to add fine tuning to the model. 
The same excuses that were used in model A can be argued here. Notice that in table I VII 
the axion couplings are identical to those of model C. This is because the splitting of the 
isospin doublets forces us to identify Q as the third generation and qi as the first and second 
generation in both models C and D. Even this model renders the FN mechanism essentially 
useless. Indeed, there are too many fine tunings so the FN mechanism is not really doing 
what we were set out to do at the beginning. Its axion falls an order of magnitude short of 
the lower sill of the axion window and one can say that this model is essentially ruled out. 

The upshot of this analysis is that for each model the axion couplings are and there 

is no accidental fine-tuning which might radically alter the lower sill of the axion window 
for each model. Thus, we find that these viable models result in one harmful, one harmless, 
and two mostly harmless axions. 



Dimension 4 couplings 


hqju + h'^Qd + hHe 


High dimension couplings 


(phQu + (j)h'^qid + {(j)hlf 


(</)) ~ 


10^ GeV 




10^ GeV 


9app 


-0.192 


9ann 


0.187 


E 


55/24 




3 X 10-1^ GeV-i 



TABLE VI: Model D 



VI. GENERAL CONSIDERATIONS AND HEAVY AXIONS 



So far, we have built four simple models of flavor textures by following the D-brane 
effective field theory paradigm. General models will be more complicated, but there are 
universal features worth discussing. First, considering the general procedure we have used 
so far, the CKM matrix is completely undetermined. This is either because the left handed 
quark doublets all appear on the same footing in the Lagrangian (the four stack models), 
or because the right handed fermions all appear on the same footing in the Lagrangian 
(the three stack models). It would be nice if the explanation for the CKM matrix was 
part of the phenomenology addressed by the Froggatt-Nielsen mechanism itself and would 
not require an independent solution. In many top-down string theory models this is done 
geometrically, the various Yukawa couplings are related to areas of triangles (see 



38 



39i for 



recent progress in the computation of Yukawa couplings in these models), but geometry by 
itself usually goes beyond effective field theory. Some of these textures might also arise from 
discrete symmetries, but we are trying to avoid that kind of argument: to our knowledge 
there is no satisfactory solution of the CKM matrix at this level. Indeed, the Frogatt-Nielsen 
mechanism is the preferred way to address this issue in phenomenology anyhow. 

The second statement that is important to understand is that by using the three and 
four stack models, independently of how the scalar sector is realized in perturbative D-brane 
models, we are committing ourselves to a low energy effective field theory with Peccei-Quinn 
symmetries that are only broken non-perturbatively. These approximate global symmetries 
are always generated by the Ibafiez-Quevedo mechanism. Thus the low energy effective field 



theory will necessarily have light axions arising from the phase factor of the order parameter 
that breaks the Froggatt-Nielsen symmetry: this can be a complicated linear combination 
of the scalars that are present, and might even belong to the Higgs sector of the low energy 
effective field theory. Many of these models are ruled out by current data and therefore 
axion constraints should be part of the analysis of top down approaches to understand the 
Standard model or its supersymmetric versions. The Ibaiiez-Quevedo mechanism works just 
as well in supersymmetric models as in non-supersymmetric models and is a property of the 
perturbative structure of open string models. 

Also notice that as far as low energy effective field theories are concerned, the important 
detail to determine the axion model is determined by the axion photon couplings (^q^ and the 
axion decay constant. Together they determine the mass of the axion. The most important 
relation for comparison to experiments like CAST is given by 

2-n \N J m^f^ 



For Grand Unified models E/N ^ 8/3 (see for example |27l|) producing a slight cancellation 
in the axion photon coupling , while for the KSVZ model [23|, l2J| E/N = 0. In our cases 
E/N — 1.92 is larger than the KSVZ model, and can be twice as large. This makes these 
models more sensitive to the constraints from experiments like CAST {34 1, although our 
considerations put the axion photon coupling for reasonable models in the range 10^^^- 
10^^^, but the precise value of the decay constant depends on values of coupling constants 
that we have set to be of order one, but as yet unspecified. If moreover axions comprise the 
dark matter halo, then our numbers can be compared with the results of ADMX jssi and 
the parameter space for our models can get constrained. 

A second point that is important to consider is that we assumed at the beginning that the 
Frogatt-Nielsen symmetry was being gauged was anomalous, and the anomaly was cancelled 
with an axion-like field that got eaten up to make a massive Z' near the string scale. This 
is a closed string field. However, the mass of the Z' receives two contributions, one from 
the Green-Schwarz anomaly cancellation mechanism, and another from the usual Higgs 
mechanism from breaking the FN-symmetry. If the mass of the Z' is dominated by the 
Higgs term, the low energy axion that survives is the one arising from the closed string 
sector. The associated axion decay constant is then bounded above by the mass of the Z'. 
This means that the closed string axion would have a relatively low decay constant and 



might fall outside the axion window. Thus in the general theory analysis it is possible to 
push down the axion decay constant, but this does not push it up. Indeed, the mass of the 
Z' fields arises from the terms 

^i^..(^;)(^^) + (0)'(^r)' (4) 

In the case the second terms dominates, the quantity K is bounded above by (0), and 
K plays the role of the axion decay constant for the closed string modes. Obviously the 
mixings complicate the analysis substantially and one also expects that the effective field 
theory might break down before the Frogatt-Nielsen scale is reached, making it necessary to 
build a more complete string model. This is outside the effective field theory paradigm we 
have been pursuing, but it is a distinct possibility. Fortunately one can not hide the axion 
this way so the models are testable. 

The axions appearing in these models can trace their origins to stringy effects, so it 
is useful to consider other stringy contributions to the axion mass. Worldsheet, D-brane, 
and other instantons can have possibly large, nonperturbative contributions to the axion 
potential. Combining the effects of stringy instantons and QCD, axions develop a potential 
of the form jlO | 

V{a) ~ -2M,^e-^-' cos(a//„ + V-) + a'^^ 

J a 

All quantities in the potential are predetermined in the models discussed here, so we can 
treat Smst as a free parameter. Ordinarily, the instanton action is bounded below in order 
for the axion to solve strong CP. This bound is^ 

, \ I 60 (model A) 

Sinst > log ' 



10 ^^mlF^J 1 125 (model B) 
In the context of the models discussed here, however, one may take the view that axions 
are simply unwanted byproducts of the model, and need not solve the strong CP problem. 
There is then freedom to consider cases where stringy axion mass contributions are impor- 
tant. If the stringy effects dominate, this divorces the relationship between mass and decay 
constant that is characteristic of the QCD axion (although the couplings all still depend 



^ If Mg is taken to be the Planck scale, Sinst ^ 200. 



on fa). If the axion can be made heavy enough, rUa > 100 MeV, this can avoid the usual 
astrophysical constraints, and save e.g. model A (models C and D should still be disfa- 
vored because of their poorly-behaved mass hierarchies) . This gives an upper bound on the 

instanton action 



For the models at hand the bounds on Sinst for QCD and heavy axions are incompatible, 
so these axions cannot simultaneously resemble QCD axions while avoiding the astrophysical 
constraints. 

However, the decays of heavy axions run the risk of interfering with big bang nucle- 
osynthesis. To avoid this, one can require that primordial axions decay away before the 
nucleosynthesis era In — 10^ — 10^ s, i.e. Ta-^xtN ^ 1- Photonic decays dominate, since 
the branching ratio to fermionic decays is 



Since Xa is of the same order as ^'077, this branching ratio is very small, peaking at 
~ 10-=^. Writing M, ~ 10"^ GeV, /„ ~ 10-^ GeV, and h ~ lO'^^ GeVs we have 



For model B, this bound matches the bound necessary to have heavy axions in the 
first place, while for model A the astrophysical bound on Sinst is more stringent, but it is 
certainly possible for these heavy axions to decay fast enough in both cases. If we think of 
these potentials as generated from an anomaly, we would need a hidden sector for which 
the U{l)pQ symmetry is anomalous, but the corresponding fermions are arranged into Dirac 
representations of the standard model. Such a model is not minimal, but it can serve as a 
field theory realization of this idea that does not require strictly stringy instantons. 





h Pa ^ Pa ^ 



^Q6m-5/-|^Q2+24g-35i„,t/2 



This places the constraint on the instanton action 




Another possibility is that heavy axions are very long-hved, and comprise dark mat- 
ter. Sikivie {sg], Is?! has argued that caustics in the dark matter halo are indicative of a 
Bose-Einstein condensate of axion-like particles. If so, then we must have Ttnow ^ 1- Un- 
fortunately for heavy axions, their large mass makes the mean lifetime is rather short. For 
axions with rria = 100 MeV, 



This problem only gets worse the heavier the axion. Trying to avoid the astrophysical 
constraint on by allowing for heavy axions can not result in heavy axion dark matter in 
the models considered here, since fa is not a free parameter. 

Model A is ordinarily ruled out if the axion is a QCD axion, but if its mass depends 
more strongly on stringy instanton effects than its decay constant, the lower bound on /„ 
can be relaxed. This is allowable for fast-decaying axions, which means that the action of 
the instantons contributing to the axion potential must be rather small, in stark contrast to 
the QCD axion. Nevertheless, it is thus possible to open the axion window wide, and allow 
model A through. 

A second problem that we have alluded to before is that the CKM matrix is not part 
of our proposal for addressing the flavor problem. We would like to comment now on how 
to remedy this situation. The sticky point in our models is that either the left or the right 
handed quarks are essentially indistinguishable from each other as far as their quantum 
numbers are concerned. Thus a model that does better needs to distinguish both the left 
handed quarks and the right handed quarks. Let us assume that we want to do this in some 
minimal way. What we need is some gauge theoretic reason to connect some left handed 
quarks preferentially to some right handed quarks and not the others. 

Again, we want to achieve this by using a gauge symmetry reason. Since the quarks are 
attached to the same U{3) node of QCD, this is the place where we can make the biggest 
difference, thus we would want an extension of the U (3) color to a bigger gauge group. This 
should be chiral, and therefore should be represented by products of U{ni) gauge groups. 
These typically break to diagonal subgroups with standard embeddings due to the D-brane 
matter content. The simplest scenario is to have a U{3) x f/(3) symmetry where the top 
and bottom quark are attached to one U{3) and the other quarks are attached to the other 




10-15 s (A) 
15 s (B) 



f/(3). We also need a scalar that breaks to the diagonal f/(3) x f/(3) — )■ U{3). The natural 
scalar is in the (3,3) representation. Notice that this introduces a new f/(l) gauge group 
that might suffer anomalies as well and similar reasonings as before (using the Ibanez- 
Quevedo argument) lead to the possibility of a second low energy axion in the dynamics. 
This extra scalar field would appear in non-renormalizable terms in the lagrangian that 
provide for the smallness of the mixings between the third family and the other two families. 
This is indeed what a deconstruction of an extra dimension for QCD would look like 40 1 
(supersymmetric extensions of the standard model that have this flavor structure have been 



considered recently in 



4l| ). This is only one step away from thinking about this problem 



using large extra dimensions as they appear in a full compactification of string theory and 
suggests that in the end this problem can always be argued to be resolved geometrically. 
Naturally, the analysis here becomes more complicated and will not be pursued further in 
this paper. 



VII. CONCLUSION 



We have considered using the Froggatt-Nielsen mechanism to generate the minimal 
amount of hierarchy between fermion masses, separating them into light and heavy tiers. 
The novelty of our construction is that the Frogatt-Nielsen symmetry has mixed anomalies 
that are cancelled with the Green-Schwarz mechanism. This anomaly cancellation provides 
a Stuckelberg mass for the Z' associated to this symmetry. The peculiarities of realiza- 
tions of this class of models in perturbative D-brane constructions leave a residue of the 
Froggatt-Nielsen symmetry at low energies that acts as a Peccei-Quinn symmetry. When 
the symmetry is broken in order to account for flavor hierarchies, a pseudo-goldstone boson 
is generated which behaves like an axion. This arises from the open string sector dynamics 
and can be considered as an open string axion. We do an effective field theory analysis of 
the physics by considering non-renormalizable contributions to the lagrangian suppressed 
by a high energy mass scale which we take to be the string scale. The ratio of the Frogatt- 
Nielsen order parameter to the string scale is determined by the quark mass ratios. The 
Froggatt-Nielsen fields can also participate in how the neutrino masses are generated, so it 
is natural to consider the neutrino masses to calculate the high scale. If the neutrino masses 
are generated by dimension five operators, the high energy string scale is of order 10^^ GeV, 



as is typical of large extra dimension models based on D-branes. If the neutrino mass is 
generated by dimension 7 operators the string scale is much lower and depends on details of 
the Froggatt-Nielsen suppressions of a given model. In general, the Froggatt-Nielsen order 
parameter serves as a benchmark for the axion decay constant and it is between one and 
three orders of magnitude below the string scale. Many possible realizations of the standard 
model within string theory could be ruled out this way, as they would lead to unacceptably 
small axion decay constants outside the allowed axion window. This issue makes it clear that 
axion bounds need to be analyzed in order to decide if D-brane models built from string the- 
ory are truly viable. Furthermore, our analysis predicts relatively large axion-photon-photon 
couplings so the models can be compared to experimental searches like those carried at the 
CAST and ADMX experiments. This gives us a prospect for ruling out such extensions of 
the standard model. 

There are ways out of this dilemma that would require additional stringy instantons to 
make the low energy axions heavier in order to avoid astrophysical constraints from nuclear 
physics processes. We also found that in the simplest models we consider we have nothing to 
say about the CKM matrix and to find an acceptable solution that includes this information 
it is likely that one needs to deconstruct extra dimensions so that strong QCD dynamics 
arises from a C/(3) x U{3) U{3) gauge symmetry breaking at high energies or an even 
bigger extension. The details of such models are beyond the scope of the present paper. 

If one is willing to give up the possibility of a string embedding, one can assign arbitrary 
U{l)a,b charges to the various fields; this allows higher-dimensional FN terms of the generic 
forms 0^ {hq^qpi) and (p"{hliy , giving masses rUq ~ {h) (^■^'^ i ""^i^ — (^) ^'^^ ^^S^^ 
quarks and neutrinos, respectively. Both M and (0) can be estimated from combining the 
information of these ratios. In particular, the most important for us is the axion decay 
constant which is essentially given by the vacuum expectation value (0) . This is given by 

{(f)) ~ mf ^-"-'^/''m^^+i^/^m-i 

where is the mass of the heaviest light quark. This is very model dependent. In general 
one obtains a lot of possibilities this way and one needs to pick the one that corresponds to 
the highest value of (0) , by assuming that the other couplings that give rise to a high scale 
have extra fine tunings that make them smaller than our naive guess would suggest. 

To conclude, we find that some string inspired minimal models for fiavor are able to 



predict axions within the allowed axion window by assuming that the scale of new physics 
(taken as the string scale) is closely related to the scale at which the neutrino masses are 
generated. The models also predict an intermediate string scale that can be rather light in 
the neutrino masses arise from high dimension operators (dimension 7). If the string scale is 
too light, the axion falls outside the allowed window. In each case, the models fail to explain 
the smallness of some of the Yukawa couplings. It would be interesting to study further how 
such models can be embedded in supersymmetric setups so that the hierarchy problem can 
also be addressed. 
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